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Abstract 



Q 1 INTRODUCTION 



We can say that quantum computation is a young and promising subject and can revolutionize our 
computational way of tackling problems, especially the simulation of quantum problems numerically. 
O" 1 ' Nowadays quantum computation is a branch of the wider subject quantum information 12 El H] 
£> . which embraces besides quantum computation, for instance, quantum communication and quantum 
cryptography. P. Benioff and R. Feynmann [Oj were the pioneers of quantum computation 
suggesting the build up of computers based on the principles of quantum mechanics. Following 
the lead of Benioff and Feynmann, D. Deutsch in 1985 [7] made concrete proposals exploring some 
properties of quantum mechanics to obtain unprecedented parallelism in computation which repre- 
sented a really breakthrough to the subject which since then has developed quickly although we can 
still say that perhaps we are far from building a working practical and economically viable hardware 
based on quantum mechanics, that is, unfortunately we are still far from having, for instance, a 
quantum PC for helping us in working out problems and, why not, for our enjoyment. However, 
we should point out that although it seems that a quantum computer would speed up the solution 
of some problems and even turn out possible the solution of some problems intractable by classical 
computation, nobody still knows what is the real power of a quantum computer when compared 
with classical computers. Despite this fact we should go on in the enterprise of developing research 
on this field which can at least give us new insights on quantum theory and quantum information 
science. 

Important quantum algorithms are the Deutsch algorithm [7J H], the Deutsch-Jozsa algorithm 
jEJII], the Simon algorithm [315, the Shor algorithm [T0| H] and the Grover algorithm [TT| Hj. 
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In this work we present an implementation of the Deutsch 0H] and the Deutsch-Jozsa algorithms 
jHlH], via cavity QED. The Deutsch algorithm was the first concrete proposal of computation making 
use of the special features of quantum mechanics. A recent alternative proposal of realization of the 
Deutsch algorithm is presented in 



2 CAVITY QED REALIZATION 

Let us start revising the Deutsch problem. Consider an arbitrary Boolean function F : {0, 1} — ► 
{0, 1}. There are four such a functions corresponding to two possible arguments and two possible 
values. For two of them F(0) = F(l) and in this case we say that F is constant. For the cases 
in which F(0) ^ F(l) we say that F is balanced. Suppose we do not know the function and we 
are given an Oracle which can evaluate it and gives us the result. Notice that in order to decide if 
F is constant or balanced we will have to use the oracle twice to know its value for and 1. The 
Deutsch algorithm can solve this problem with just one call of the oracle. Let us see how the Deutsch 
algorithm works. First let us suppose we have a gate (F— gate). The action of the F— gate is 



\x,y) — >\x,y® F(x)). 



(2.1) 



The Deutsch algorithm employs a F— gate with F being our function which we want to decide if is 
constant or balanced. Considering that the input two qubit state is 



the output is 



l^n) = 1^)^(1°) " I 1 )) = \ (M> - \X, 1)), 



|VW> = -j={\x, F(x)) - \x, 1 © F(x))). 



Since f{x) = or 1 this can be written as 



|Vw> = (-i) F( ^(M> - |s, 1)) = (-l) F ^\x)^=(\0) - 



(2.2) 



(2.3) 



(2.4) 



and as we see the effect is to change the state of the \x) qubit to (— 1) F( - X ^\x). Therefore the value of 
the function is in the phase of the state \x). 

The actual circuit of the Deutsch algorithm is shown in Fig. 1 where if is a Hadamard gate 



H — — = 



1 

V2 



1 1 
1 -1 



(2.5) 



and 



|0> 
|1> 



V2 



(|0> + |1», 



^(10) -ID). 



Then, as we see, after the action of the F— gate we have 

io, i) — ^[(-ir (o) io> + (-if (i) ii>]-^(io> 



(2.6) 
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and, for F(0) = F(l), the upper qubit is therefore in the state ±4=(|0) + |1)) and after the application 
of the last Hadamard gate it will be in the state ±|0). For F(0) ^ F(l), the upper qubit is therefore 
in the state ±^(|0) — |1)) and after the application of the last Hadamard gate it will be in the 
state ±|1). Thus a single call to the quantum oracle followed by a measure of the upper qubit in the 
computational basis gives us the answer of the problem. Notice that the state in which the lower 
qubit is left is not important. 

Let us see now how we can implement experimentally the Deutsch algorithm. We start assuming 
that we have a cavity C prepared in the state 

(|Q) C -|1) C ) 

\-)c = ^ • (2.T) 

In order to prepare this state, we send a two-level atom AO, with |/ ) and |eo) being the lower and 
upper level respectively, in the state 

I V>>ao = I eo)+ I /o», (2-8) 

through C, for ^40 resonant with the cavity. If g is the coupling constant and r the atom-field 
interaction time, under the Jaynes-Cummings dynamics for gr = 7r/2, we know that the state 
|/o)|0)c does not evolve, however, the state |eo)|0)c evolves to —i\fo)\l)c- Then, for the cavity 
initially in the vacuum state \0)c, we have 

m^ l0)c ^ M m^m =lfo)Hc (2 . 9) 

Then let us assume that we have a cavity C prepared initially in a state | —)c and an atom Al 
prepared initially the state 

I = ^(1 I (2.10) 

Consider now that for an atom Ak, with \f k ) and \e k ) being the lower and upper level respectively, 
the transition | ff.) ^| e k ) is far from resonance with the cavity central frequency as shown in Fig. 
2. In this case the effective Hamiltonian is given by [Hj 

2 2 

H = h 9 -{a)a + 1)(| e k )(e k | -h^ata \ f k )(f k |), (2.11) 
and the time evolution operator is given by JH] 

U(t) = e -^ a+ V | ek )( ek | +e ^« | |, (2.12) 

where (p = g 2 r / A, g is the coupling constant, A = u e — cjf — u> is the detuning where u e and Uf are 
the frequencies of the upper and lower levels respectively and uj is the cavity field frequency and r 
is the atom-field interaction time. Now, we are going to send Al through cavity C where the atom 
interacts dispersively with C according to ([2. 12)1 . Let us take tp — 7T. Then, after Al flies through C 
prepared in state (|2.7jh we have 



1 

)ai-c = 2' 



|ei> + |/i»(|0>o+|l>c). (2-13) 
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If (p — 27T we have 

\1>)AL-c = \{\ei) + \fi))(\0)c-\l)c). (2.14) 

Now, if we use the notation 



ejt) =| 0)aa:, 

I /*)=|1)a*, (2.15) 

we can rewrite ()2.13|) and (j2.14j) as 

W)ax-o = ±^[(-l) m \0) A1 + (-l)W\l) M ]±=[\0) c - (_i)^o)^(i)| 1)d (2 . 16) 

Now we make use of the Hadamard gate H ()2.5|) . Then, in the case of F(0) = F(l) the atom will be 
in the state ±A^(\0)ai + |1)ai) or ±^=(|ei) + |/i)) and after we apply the H gate we get ±|0)ai or 
|ei). In the case of /(0) ^ /(l) the atom will be in the state ±4=(|0).ai — |1)ai) or ±-^(|ei) — 
and after we apply the H gate we get ±|1)ai or \f\). Notice that the state ^[|0) - ( — l)^(°)®^Ci) 1 1> ] 
in which the cavity is left is not important. 

Another possible implementation of the Deutsch algorithm is possible considering a three-level 
cascade atom Ak with | e k ), | /&) and | g^) being the upper, intermediate and lower atomic states. 
As above, we assume that the transition | /&} ^| e&) is far enough from resonance with the cavity 
central frequency such that only virtual transitions occur between these states (only these states 
interact with field in cavity C). In addition we assume that the transition | e k ) ^| gj.) is highly 
detuned from the cavity frequency so that there will be no coupling with the cavity field (see Fig. 
3). Here we are going to consider the effect of the atom- field interaction taking into account only 
levels | fk) and | g k ). We do not consider level | e k ) since it will not play any role in our scheme. 
Therefore, we have effectively a two-level system involving states | fk) and \gu). Considering levels 
| fk) and | gk), we can write an effective time evolution operator (see ([2.12)) ). 

U k (t)=e^ a \f k )(f k \+\g k )(g k \. (2.17) 

A coherent state \a) is obtained applying the displacement operator D(a) = e (aat_Q * a ) to the vacuum, 
that is, | a) = D(a)\0), and is given by 

\a) = e-^Y.^k\n)- (2.18) 
n=0 v n\ 

We define the even and odd coherent states as 

1+) = -^(|a> + |-a» = -^(|0> c +|l>a), 

I-) = ^=(|a>-|-a)) = -^(|0>o-|l)a), (2-19) 

with N ± = 2 (l ± e _2 ' a ' 2 ^ = 2 and (+ | — ) = 0. Now, let us see how we can prepare a state |— ). 
Suppose we prepare cavity C initially in a coherent state | — a). Then we prepare a two-level atom B, 
with | /) and | g) being the upper and lower state respectively, in a coherent superposition, sending 
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B in the lower state 
according to 



and we get 



through a first Ramsey cavity Rl where the atomic states are rotated 



C 9 C f 
~ C f C 9 



(2.20) 



I $)b = Cf\f) + c g \ g). (2.21) 

After that, B flies through cavity C and, taking into account the time evolution operator (|2.17jl . 
after B pass through C the state of the system B — C, for ip = ir, is given by 

| ip) B -c = Cf\f)\a) + c g \g)\ -a) 

Then, we send B through a second Ramsey zone R2 where the atomic states are rotated according 
to 



R2 = ^= 



1 

V2 



1 -i 
-i 1 



(2.22) 



that is, 



/> 
9) 



1 

72 
1 

71 



(I />-*!</», 
H\f)+\g)), 



(2.23) 



and therefore, the state of the system B — C will be 

\*I>)b-c = ^[( c / - ic a) I +) + ( c / + ic 9 ) I -)] I /) 

+ ^[H C / + c g) I +)) - ( ic f + c g) I -)] I 9), 

Now, in order to obtain a state |— ) in cavity C, we detect atom B in | /) or in | g). If we detect | /) 
we must have have c/ = ic g and if we detect | g) we we must have have c/ = — ic g . 

Then, we start assuming that the cavity is prepared in a state |— ) and atom Al is prepared in a 
state 



Al 



1 

V2 



(IA>+I<7i»- 



After the atom fly through the cavity, for (f = n, we have 

1 



If ip = 27T we have 



Now, if we use the notation 



Ai-c = ^-\h) + \gi))(\o)c-\i)c). 



)Al-C=^(\fl) + \9l))(\0)c-\l)c). 



(2.24) 



(2.25) 



(2.26) 



fk) — I Q/Aki 
9 k) =| 1)a&) 



(2.27) 
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we can rewrite (|2.25j) and (|2.26j) as 

m A i-c = ±^[(-l) F(0) |0) A1 + (-l)^)|l) A1 ]i=[|0) c - (2.28) 

Then, in the case of F(0) = F(l) the atom will be in the state ±4^(\0} A1 + \l) A i) or ±^(|/i) + \gi)) 
and after we apply the H gate we get ±|0)ai or \fi). In the case of F(0) ^ F(l) the atom will be 
in the state ±-^(\0)ai — 1 1) ai) or ±^(|/i) — \gi)) and after we apply the H gate we get ±|1)ai or 
\9i)- 

The speed up achieved in the Deutsch algorithm over the classical algorithm is just a factor of 
two. However, in a generalized problem, we are going to discuss now, the speedup is much more 
than in the classical case showing the power of the quantum computation. Let us examine now the 
Deutsch- Jozsa algorithm to solve the generalized Deutsch problem. Let F : {0, 1}®" — > {0, 1} be 
a Boolean function of a n-bit integer argument and we assume allow only those / that are either 
constant or yield for exactly half of the arguments and 1 for the rest, that is, it is balanced. Given 
an oracle that evaluates the function for a given argument, the problem is to decide if F is constant 
or balanced. There are 2 n possible arguments and to solve the problem classically we will have to 
calculate the function for 2 n_1 + 1 arguments in the worst case. Then, we see that the computational 
resources required to solve the problem grows exponentially with the bit size n of the input. The 
Deutsch- Jozsa algorithm can however solve this problem very easily. This algorithm makes use of a 
quantum /—gate that is a generalization of the one we used above in the Deutsch algorithm. The 
action in the computational basis is similar to that of the ordinary F— gate, except that \X) here 
is a computational basis state of a n— bit register. If the bottom qubit is in the state ^-(|0) - |1)), 
then the state of the upper register is transformed according to 

\X) — > (-1) F W\X). (2.29) 

The quantum circuit of the Deutsch- Jozsa algorithm is shown in Fig. 4. The upper input is |0)|0)...|0) 
(with n factors) and H® n — H ® H eg) ... ® H. Then, the effect of the first Hadamard gates on the 
top input state is 

HnO>|0>...|0> = -/L^ =o V>. (2.30) 

The Hadamard gate acting on the bottom input state |1) produces ^(|0) — |1)). Thus the /—gate 
changes the state of the register to 

ii:i::(-i) F,x, iA->. (2.3D 



Now notice that the action of the Hadamard gates on a computational basis state is given by 

1 ■*-^2 n -l 



H*»\X) = —YZ^(-l) XJr \Y). (2-32) 



where X.Y is the bitwise scalar product where for X = x n _i...xix and Y = y n -i---VxyQ we have 
X.Y = ®2^o x iVi an d © is addition mod 2. Then, the state of the register will be 

1 J2to(-^ f(X) H^\X) = ^E^E^-l^^in (2-33) 



Now, the amplitude of \Y — 0) = |0)|0)...|0) is ^kY.x=o(~ and if F is constant this gives us 
±1. On the other hand, if / is balanced one half of the terms in the sum cancels exactly the other 
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half and the result is 0. Then, the probability of observing \Y — 0} = |0)|0)...|0) is 1 if F is constant 
and is if it is balanced. 

Let us see now how we can implement the Deutsch-Jozsa algorithm experimentally. Considering 
atoms Al,A2,..., An prepared in a state like ()2.10|) . After they fly through C prepared in state (|2.7p . 
taking into account (J2.12)) for (p = n we have 

W)ai-c = ^(-h> + (-lea) + l/2»-^Hen> + l/»»^(|0) O - (-l) n |l)c). (2.34) 

If ip = 2tc we have 

WiAi-c = ^(|d> + l/i))^(|e 2 ) + |/ 2 »...^(|e»> + |/»»^(|0> o - \l)c). (2.35) 



If we assume atoms Al, A2, An prepared in a state like (|2.24J) and cavity C prepared in state 
|— ), taking into account 1)2.17)1 for <p = n we have 

W)AX-G = ^f(-IA) + l^l»^(-i/2> + \92})...^Hfn) + k»^(|0>C " |1> C ). (2-36) 

If ip = 2tt we have 

\1>)ai-c = ^(l/i) + l5i))^(l/ 2 > + |a,»...^(|/„> + |»»»^(|0>o " |l>o). (2-37) 
Now, if we use the notation JZH51) and $TT7\ we can rewrite (J23E1), (jZHBD and (pOTjl as 

1 ^2 n -l 



>A1-C 



±^ f e;: (-i) f(x) i^)ai^)c (2.38) 



where, in the case of (j2.34j) and ()2.35|) . if F is balanced, 

l^)c = ^ f (|0)c-(-ir|l)c), (2.39) 

and if / is constant, 

|^)c=^(|0)c-|l)c). (2.40) 

In the case of fl2~3fijl and (f2~5?j) . 

I^)c = ^(|a> - | - a)) = ^(|0) c - |l>o). (2.41) 
Now notice that the action of the H gates on a computational basis state is given by 

on i 

H® n \X) A = "^EyJo (-l) Xy |>V (2.42) 
Then, the state of the register will be 

1 E^C-i^^Wa = *YZ^YZ^(-i) fw+x ' y \y)* ( 2 - 43 ) 



Now, the amplitude of |0)^ = |0}^i|0)^2...|0}^ n is i^Yyx=o(~ l) F< - x - ) and if F is constant this gives 
us ±1. On the other hand, if / is balanced one half of the terms in the sum cancels exactly the 
other half and the result is 0. Then, the probability of observing \0)a = |0)xi|0)yi2"-|0)An is 1 if F 
is constant and is if it is balanced. Therefore, if we detect |0)^i|0)^ 2 ---|0)^ n then F is constant. 
Concluding, a single call to the quantum oracle followed by a measure of the register and checking 
if the result is |0)Ai|0) J 42---|0)An allows us to decide if the function is constant or balanced and we 
have achieved an exponential speedup over the classical computation. Finally, let us analyze the 
feasibility of the experimental implementation of the above algorithms. Considering Rydberg atoms 
of principal quantum numbers 50 or 51, the radiative time is of the order of 10 -2 s and the coupling 
constant g is of the order of 2n x 25 kHz ^SICSIEI and the detuning. A is of the order of 2ir x 100 
kHz. Taking into account that <p = g 2 r/A, for ip = it we have an interaction time r = 8 x 10~ 5 s 
and we could, in principle, assume a time of the order of 10~ 4 s to realize the algorithm which is 
much shorter than the radiative time. We have to consider also the cavity decay time which in recent 
experiments, with niobium superconducting cavities at very low temperature and quality factors in 
the 10 9 — 10 10 range, have a cavity energy damping time of the order of 10 to 100 ms, and which 
could be larger than the required time to perform the algorithm. 
Figure Captions 

Fig. 1 - Quantum circuit of the Deutsch algorithm. 

Fig. 2 - Energy states scheme of a two- level atom where |e) is the upper state with atomic 
frequency u e , \f) is the intermediate state with atomic frequency uif, and u is the cavity field 
frequency and A = {uj e — ujf) — uj is the detuning. The transition | /) ^| e) is far enough of 
resonance with the cavity central frequency such that only virtual transitions occur between these 
levels. 

Fig. 3 - Energy states scheme of a three-level atom where |e) is the upper state with atomic 
frequency u e , \f) is the intermediate state with atomic frequency u>f, \g) is the lower state with 
atomic frequency uj g and uj is the cavity field frequency and A = (u e — ujf) — uj is the detuning. The 
transition | /) ^| e) is far enough of resonance with the cavity central frequency such that only 
virtual transitions occur between these levels (only these states interact with field in cavity C). In 
addition we assume that the transition | e) g) is highly detuned from the cavity frequency so that 
there will be no coupling with the cavity field in C. 

Fig. 4 - Quantum circuit of the Deutsch- Jozsa algorithm. 
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